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Abstract 

We examine how the production of TeV scale black holes at colliders is influenced 
by the presence of Lovelock higher-curvature terms in the action of models with large 
extra dimensions. Such terms are expected to arise on rather general grounds, e.g., 
from string theory and are often used in the literature to model modifications to the 
Einstein-Hilbert action arising from quantum and/or stringy corrections. While adding 
the invariant which is quadratic in the curvature leads to quantitative modifications in 
black hole properties, cubic and higher invariants are found to produce significant qual¬ 
itative changes, e.g., classically stable black holes. We use these higher-order curvature 
terms to construct a toy model of the black hole production cross section threshold. 
For reasonable parameter values we demonstrate that detailed measurements of the 
properties of black holes at future colliders will be highly sensitive to the presence of 
the Lovelock higher-order curvature terms. 
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1 Introduction 


The introduction of large extra dimensions by Arkani-Hamed, Dimopoulos and Dvali^] offers 
us the possibility that the true fundamental scale of gravity, M*, may not be too far above 
the weak scale, ~ TeV. In this scenario, gravity is allowed to propagate in all D = 4 + n 
dimensions while the Standard Model(SM) fields are conhned to live on a three-dimensional 
’brane’ which we assume to be flat. One then hnds M* is related to the usual 4-d (reduced) 
Planck scale, Mpi, via the now famous expression 

m% = v,m:*\ ( 1 ) 

where Vn is the volume of the compactihed extra dimensions. These compactihed dimen¬ 
sions are usually assumed to be flat, i.e., they form an n-dimensional torus so that, if all 
compactihcation radii, Rc, were the same, Vn = (27ri?c)”- This simple ADD picture has 
three basic predictions that have gotten significant attention in the literature [ 2 ] over the 
past few years: {i) the emission of graviton Kaluza-Klein states during the collision of SM 
particles leading to hnal states with apparent missing energy j2l 0 Ej; (ii) the exchange of 
graviton Kaluza-Klein excitations between SM helds leading to new dimension-8 contact 
interaction-like operators with distinctive spin-2 properties[21 0 0; {iR) the production of 
black holes(BH) at colliders and in cosmic rays with geometric subprocess cross sections 
once energies greater than ~ M* are exceeded [3 0 0 [10] . While (i) and (ii) result from an 
expansion of the gravitational action, i.e., the D—dimensional Einstein-Hilbert(EH) action, 
to leading order in the gravitational field and are in some sense perturbative, {Hi) on the 
other-hand relies upon the full non-perturbative content of the EH action. Here one is really 
testing D—dimensional General Relativity and not just the ADD picture. Within the ADD 
scenario, collider measurements of (i) and {ii) type processes will tell us the values of n and 
M*|2|. 
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If anything like the ADD scenario is realized in Natnre it will have to be part of a ranch 
larger framework, i.e., ADD is at best an effective theory that operates at energies below the 
scale M*. It is reasonable to expect that at least some aspects of this more complete theory 
may leak down into the collider tests listed above and may lead to signihcant qnantitative 
and/or qnalitative modihcations that can be probed experimentally. Here we examine one 
snch extension of the basic model: the existence of higher cnrvatnre terms which angment the 
EH action. Snch terms are expected to be present at some level on rather general gronnds 
from string theory [TTl IT^ or other possible high-scale completions of General Relativity. 
These terms can also be thonght of as being generated by qnantnm corrections to the ordinary 
EH action. Given arbitrary powers and derivatives of the cnrvatnre tensor there are a hnge 
nnmber of possible higher-order invariants from which to choose. In order to decide what to 
add to the EH action we need some form of gnidance in making possible selections. 

Fortnnately, a certain special class of snch invariants with very interesting proper¬ 
ties was hrst generally described long ago by Lovelock^H] and, hence, are termed Lovelock 
invariants. These are constrncted ont of powers of the cnrvatnre tensor with no additional 
derivatives. They are also sometimes referred to in the literatnre as generalized Enler den¬ 
sities since their volnme integrals are related to the Enler characteristics. These Lovelock 
invariants themselves come in hxed order, m, which we denote here as that describes 
the nnmber of powers of the cnrvatnre tensor, contracted in varions ways, ont of which they 
are constrncted. Apart from normalization factors we can express the Cm as 


Cr 


m p CiDi D 

.nt 


AiSi 


^AmB- 


Cm Dm 


m^m 


( 2 ) 


where ^clEl cZdZ totally antisymmetric prodnct of Kronecker deltas and Rab 

is the D-dimensional cnrvatnre tensor. Fortnnately, as can be seen by this dehnition, the 
nnmber of snch invariants that can exist in any given dimension is highly constrained. For a 
space with an even nnmber of dimensions, D = 2m, the Lovelock invariant is a topological 
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one and leads to a total derivative, i.e., a snrface term, in the action. All of the higher order 
invariants, D < 2m — 1, can then be shown to vanish identically by various curvature tensor 
index symmetry properties. On the otherhand, for the cases with D > 2m + 1, the Cm are 
true dynamical objects that once added the action can signihcantly alter the held equations 
normally associated with the EH term. However it can be shown that the addition of any 
or all of the Cm to the EH action still results in a theory with only second order equations 
of motion as is the case for ordinary Einstein gravity. Furthermore, variation of the new 
action leads to modihcations of Einstein’s equations by the addition of new terms which are 
second-rank symmetric tensors with vanishing covariant derivatives and which depend only 
on the metric and its hrst and second derivatives, he., they have the same general properties 
as the Einstein tensor itself but are higher order in the curvature. These properties are quite 
special. Generally, the addition to the action of arbitrary invariants formed from ever higher 
powers of the curvature tensor will lead to equations of motion of ever higher order, i.e., 
ever more co-ordinate derivatives of the metric tensor and graviton held, e.g., terms with 
quartic derivatives. Such theories will have very serious problems with both the presence of 
ghosts as well as with unitarity^J. The Lovelock invariants are constructed in such a way 
as to be free of these problems making them quite exceptional. Interestingly the Lovelock 
invariants are found to be just the forms taken by the higher order curvature terms generated 
in perturbative critical string theory HU Hg. This is perhaps what we might have expected 
if string theories are to avoid these troublesome ghost and unitarity issues. 

In the case of the ADD model, we are generally considering the possibility that the 
number of extra dimensions lies in the range 2 < n < 6.^ We now imagine extending the D- 
dimensional EH action of ADD to include all of the allowable Lovelock invariants for a given 

value of n with arbitrary coefficients. What do these additional terms do? Firstly, as we 

^Note that the case n = 1 is excluded by Solar System measurements and the value n = 2 with a low M* 
is somewhat disfavored by astrophysical constraintsP). 
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will see below based on simple dimensional analysis, these new invariants are suppressed in 
comparison to the usual EH action by powers of and so are only important as TeV-scale 
energies are approached. Secondly, a short amount of algebra demonstrates that their impact 
on the ADD model itself is rather subtle since the compactihed space is assumed to be flat 
and the SM helds are restricted to a 3-brane. In particular the leading-order expression for 
the interaction of the graviton Kaluza-Klein modes, with the localized SM helds remains 
as in the EH case: 

S.„, = / dV cPy (j/), (3) 

where is the stress-energy tensor for the SM helds and with y labeling the coordinates of 
the compactihed torus. (Graviton self-interactions, not generally probed in the ADD scenario 
are modihed by the Lovelock terms.) The mass spectrum of the graviton KK states is also 
unchanged. This means that the signals (i) and {ii) remain unaltered by these modihcations 
to the ADD action. Note that if the geometry of space was curved, as in the Randall- 
Sundrum model m, the couplings of graviton Kaluza-Klein excitations to SM matter would 
be altered jini by the presence of the Cm as would the KK spectrum. It is interesting to note 
that in the RS case we observe that as the curvature parameter k/Mpi —>■ 0 we hnd that the 
matter-graviton couplings return to those of usual RS scenario even though £2 is present 
with a non-zero coefficient. What about {Hi), the production of BH at colliders? This is the 
subject of the current paper as presented in our discussion below. 

The outline of the paper is as follows: In Section 2 we will present the basic back¬ 
ground and formalism associated with the Lovelock invariants and show how they would alter 
EH expectations for the Schwarzschild radius, temperature, entropy, free energy and heat 
capacity/specihc heat of BH. In Section 3 we will make a numerical analysis of the Lovelock 
induced modihcations for TeV scale BH and determine how their anticipated properties at 
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colliders would be affected. In particular we will show that while the presence of £2 can 
quantitatively alter the details of the properties of BH, the presence of £ 3^4 terms may lead 
to signihcant qualitative modihcations in their properties, e.g., the existence of classically 
stable BH with minimum allowed masses and corresponding mass thresholds. In Section 4, 
based on the analysis in the previous sections, we present a toy model for the threshold be¬ 
havior of the BH production cross section at colliders in odd dimensions. Section 5 contains 
a discussion and our conclusions. 


2 Formalism 

Before beginning the general discussion, the hrst question one might ask is ‘Just how many 
Lovelock invariants are there for our cases of interest?’ For simplicity, let us begin with 4-d; 
apart from numerical factors, only £0, which we can think of as the cosmological constant, 
and £1 = R, the usual Ricci scalar of EH, can be present and dynamical. The invariant of 
the next order, £ 2 , can be identihed with the Gauss-Bonnet(G-B) invariant, —AR^^R^^R 
RfivagR^^"^^ 1 which is a topological term as well as a total derivative and can be written in 
4-d as ^^Rpa ■ All higher order invariants can be shown to vanish. Thus the 

most general Lovelock theory in 4-d is just EH plus a cosmological constant which is just 
ordinary General Relativity. Now in 5-d, all of the £m >3 still vanish as in 4-d but the G-B 
invariant is no longer a total derivative and its presence will modify the results obtained 
from Einstein gravity. The generalization is now quite clear: for D = 5,6 only £ 0-2 can be 
present. Eor D = 7,8 only £0-3 can be present while for = 9,10 only £0-4. Since ADD 
assumes that the compactihed space is flat the coefficient of £0 is taken to be zero and, to 
reproduce conventional results, £1 is normalized so that it can be identihed with the usual 
EH term. Thus there are at most three new pieces to add to the EH action and so the 
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general form for the extended ADD model we consider is given by 


S = J d^+^x ^ 



( 4 ) 


where a, (3 and 7 are dimensionless coefficients which we take to be positive in our discussion 


below. (The size of these coefficients will be discussed later.) It is important to note that 
the fundamental mass parameter, M* appearing in the above action is the same as the one 
appearing in the ADD M* — Mpi relationship Eq.(l) and in the coupling of the SM brane 
helds to the graviton excitations Eq.(2). 

Note that M* can be related to several other mass parameters used in the liter¬ 
ature. The Planck scale employed by Dimopoulos and LandsbergjS] is given by Mpp = 
( 87 r)^/("-+^)M* while that of Giddings and ThomasjH] is found to be Mgt = 
moreover, Giudice et aL[2j employ a different scale Mp = M^. Further note that 

M* is thus correspondingly smaller than all of these other parameters with consequently far 
weaker bounds[2]. For example, if n = 2(6) and Mp = 1.5 TeV then M* = 0.60(0.38) TeV; 
collider bounds on M* are thus well below 1 TeV at present. The advantage of employing M* 
is that it a common parameter occurring not only in the original action but in the resulting 
SM matter couplings to gravity as well as in the important ADD relationship Eq.(l). 

As discussed above, in D = 4 -|- n dimensions £2 is given by the generalized Gauss- 
Bonnet form: 


r _ p2 A p I t? t)ABCD 

L,2 — it — AUABii + JXABCDJX 


(5) 


while for £3, hrst given explicitly by Mhller-Hoisson [TT)] . one has 


L,3 — K — iznn Kab + KcdEFJA AB + ort CEJA BEJA ad 


_l_ c\AT}ABCD-n t}E i Q P p^l^CT) p _L 9/1 F? F? 

-h Z4Ft KcdBeH a + oJtii itCDAB + KbaUdC 


+ 16 R^^RbcR^'A- 


( 6 ) 
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£ 4 , which has 25 terms, was hrst given explicitly by Briggs^Tj and will not be reproduced 
here but the general qualitative pattern is now obvious. As discussed above the Lovelock 
invariants lead to an augmented set of Einstein’s equations of the form 


Rab 


:9ab 


2a 


'( 2 ) 

MB 


+ .... = 


Mn+2 


T, 


AB 


(7) 


with, e.g., given by 

Rab = RRab — RacghR^^^ b — ‘^RagbrR^^ — ‘2RagR^ b + ■^9abR2 ■ (8) 

As mentioned above the set of tensors are symmetric and have vanishing covariant 
derivatives. 

The possibility that TeV scale black holes(BH) may be a copious signal for extra 
dimensions at future colliders has been discussed by a number of authors|Hl E]; the corre¬ 
sponding production by cosmic rays has been considered as well, e.g., in Ref.|18j. A hrst 
step in any consideration of this process is to ask just how large the cross section can be. A 
leading approximation for the subprocess cross-section for the production of a Schwarzschild 
BH of mass Mbh is simply given by the geometric BH size |19[0nilTnj once we reach energies 
in excess of M*, i.e., 

a ~ e{^fs - Ml ), (9) 

with s = and where from this point forward R is the (4 -|- n)-dimensional Schwarzschild 
radius corresponding to the mass Mbh- Notice that the rather unphysical threshold behavior 
of this cross section is just a simple step function. While an understanding of this kinematic 
threshold regime is likely beyond any semi-classical approach we can be sure that the true 
threshold will not be so trivial. It is now believed that this cross section expression correctly 
describes the BH production process for Mbh » Tf* up to an overall factor of order 
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unityjTH] within D-dimensional General Relativity based on the EH action. In fact, using 
this expression the production of BH at the LHC has been studied in detail, including some 
detector effects, in Refs. EH 123 . These authors have shown how it is possible to extract 
information about the values of Mbhi n, and M* from LHC data, once BH production is 
observed, owing to the very large statistics that is expected. 

The usual description of the BH production process assumes that the initial horizon 
forms around a BH state which is in general both charged and rapidly rotating due to the ini¬ 
tial large impact parameter involved in the collision and the nature of the incoming partons. 
The charge and angular momentum are rapidly shed during the ‘balding’ phase leaving a 
Schwarzschild BH which then decays by Hawking radiation described by the Hawking tem¬ 
perature, Th. (It is assumed that passing through the balding phase does not alter the 
applicability of the cross section formula above by more than factors of order unity.) Other 
global thermodynamic quantities such as the entropy. S', specific heat/heat capacity, C, and 
(scaled) free energy, F = {Mbh —ThS)/M^,, are also useful in describing the detailed nature 
of the BH. In the EH case these quantities depend solely on n and Mbh and have rather 
simple behaviors. For completeness we note that if only the EH term is present one hnds 
the following ‘standard’ relationships between these quantities which are simple monotonic 
functions of R or Mbh'- 
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( 10 ) 





where m = Mbh/M^, xq = and the constant c is given by 


{n + 2 ) 7 r(’^+ 3)/2 

f(S±3) 


(11) 


Note that all the qnantities on the left have been dehned to be dimensionless throngh ap¬ 
propriate scalings by M*. Also note that the constant c differs here from that given by 
other anthors in the literatnrejHl IHl dne to onr nse of M* as the common parameter. The 
index ‘ 0 ’ for the qnantities in this eqnation labels the fact that they arise in conventional 
4 -|- n-dimensional General Relativity assnming an asymptotically flat space and in the limit 
that corrections arising from the hnite size of the extra dimensions, of order (R/Rc)'^, can 
be neglected. Rq is thns the Schwarzschild radins for a BH of mass Mbh assnming the EH 
action m D = A + n dimensions. 

One can easily imagine that if any of the higher Lovelock terms are present then the 
conventional mass-radins-thermodynamical relations above can be signihcantly altered. Of 
conrse, hrst one needs to show that the asymptotically flat, D-dimensional Schwarzschild 
solntion still exists when the new Lovelock terms are present in the action; fortnnately this 
issne was dealt with long ago [23] and snch solntions are now well known. Since that time there 
have been many discnssions in the literatnre abont the properties of BH in Lovelock-extended 
gravityj23] hnt very little attention has been given to how snch terms may inflnence TeV- 
scale BH and the prodnction of snch BH at colliders|2ni- After some work the analyses given 
in Refs. cniEn allow ns to extract the complete expressions for the varions thermodynamic 
qnantities above for TeV-scale BH based on onr generalized action. Onr goal is to be able 
to inpnt the BH mass, he., m, as well as the valne of n and the Lovelock parameters(a, (5, 7 ) 
and then extract valnes for all of the other qnantities of interest. 

The most important relation is the one between Mbh and the Schwarzschild radins. 
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R in this more general case; one obtains: 


m = c 


+ an{n + l)x"' ^ + f3n{n + l)(n — 1)(?7, — 2)x 


n—3 


+ -yn{n + 1)(?7, — 1)(?7, — 2){n — 3)(n — 4)a; 


n—5 


( 12 ) 


where x = M^R and the ‘0’ index has now been dropped. Again, we remind the reader that 
the term proportional to 7(/9, a) is only present if n > 5(3,1). Of course what we really want 
to know is x{m) so that we must find the roots of this polynomial equation; this must be 
done numerically in general except for some special cases. Fortunately, for the parameters 
of interest to us below, we find that this polynomial has only one distinct real positive root. 
The roots of this equation do lead to some remarkable properties (that will be discussed in 
full detail below). To whet our appetite, consider a simple situation with n = 3 but where 
a = 7 = 0; in this case, the equation above can be inverted trivially. Then, for physical 
values of x it is clear that m must be bounded from below, i.e., unless the mass is in excess of 
a certain minimum, no horizon will form. This leads to a remarkable threshold-like behavior 
in the BH cross section that we will return to below. 

Once x{m) is known, the BH temperature can be obtained as usual from the metric. 
The asymptotically flat u-dimensional BH metric in the presence of the Lovelock terms has 
the general form|2S| 

^ - e~‘^^dr^ - r'^dVt^ , (13) 

with —>■ 1 as r —> oo. The dimensionless BH temperature is then given bv j23l 1^ 

r = ^(e-7'ls. (14) 

with the prime here denoting differentiation with respect to ‘r’ and the index implying 
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evaluation at the Schwarzschild radius, so that one obtains 


{n + 1) U{x) 
dvr V (x) ’ 


where 


[/(x) = x^-h an(n — l)x^ +/3n(n — l)(n — 2)(n — 3)x^ 

+ 7 n(n — l)(n — 2)(n — 3)(n — 4)(?7, — 5) 

x^ + 2an{n + l)a;^ + 3f3n{n + l)(n — l)(n — 2)x^ 


V{x) = 


X 


+ 47n(n + l)(n — l)(n — 2)(n — 3)(?7, — 4) 


( 15 ) 


(16) 


The BH entropy can then be calculated using the familiar thermodynamical relation 





dx T-^ 


dm 

dx 


(17) 


which yields 


^ = 


Attc r 


^n +2 2a[n + l)(n + 2)a;"' + 3(3n{n + l)(n + 2)(n — l)x 


.n—2 


n + 2 L 

+ 4777 ,( 77 , + 1)(77 + 2) (77 — 1)(77 — 2) (77 — 3)a;"“^ 


(18) 


Here we have followed the work in Ref. [21] and required the entropy to vanish for a zero 
horizon size. Using the expressions above the (scaled) free energy can now be calculated via 
the relation F = m — TS while the heat capacity/specihc heat is given by 


C = 


dm dm 
dT dx 


dT 

dx 


-1 


= -Attcx^-^ V 


V'lV - u'/u 


-1 


(19) 
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Figure 1; Influence of a 7 ^ 0 on the BH mass-Schwarzschild radius relationship for n = 2,4, 6 
and 20 , corresponding to the red, green, blue and magenta sets of curves, respectively. 
The solid, dashed and dotted curves in each case correspond to m = Mbh/M^ = 2,5,8, 
respectively. Here /? = 7 = 0. 

with U, V given above and with the prime now denoting differentiation with respect to x. 

Knowing the Lovelock generalizations of all of these quantities we can now proceed 
with our analysis. We will see that signihcant departures from the expectations of the EH 
action are easily obtainable. While examining these various quantities it is important to recall 
that which BH properties are measurable within the ADD framework. As mentioned above, 
collider measurements of other ADD signatures will tell us the values of n and j2j, while the 
BH mass itself can be determined by direct reconstruction [211122]; similarly the cross section, 
he., R can be obtained directly. Given Mbh-,^ and M*, various BH decay distributions 
will help to determine other quantities such as the temperature [211122] • Measurement of 
the average multiplicity in the BH decay can give us a handle on the BH temperature as 
< N >~ m/2TjS| (before fragmentation effects are included). Of course, quantities such as 
the entropy and free energy are not directly measurable but are of theoretical interest. The 
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specific heat of the BH can be inferred in some cases as it is indirectly probed by the BH 
lifetime. As we will see below, in contrast to the usual EH scenario, the BH may now have a 
positive specihc heat making it grow cooler as it sheds mass via Hawking radiation. In such 
a case the BH may have a long lifetime. 

3 Analysis 

Given the generalized expressions above we can now ask how much of a numerical influence 
the Lovelock invariants will have on the quantities of interest. The hrst issue to address is 
the values of the parameters a, (3 and 7 . As a hrst guess one might naively expect these to be 
of order unity. However, if we think of the Lovelock terms as arising from a perturbative-like 
expansion of the full action, as in string theory, the coefficients must grow smaller for the 
higher order Lovelock terms. Also, when we examine the expressions above we see that for 
any kind of perturbative expansion to make sense we must have, at least approximately, 
and yn® < 1. Since n can be as large as 6 within the usual ADD scenario, very 
crudely, we might expect that a ~ 10“^, (3 ~ 10“^ — 10“^ and 7 ~ 10“® with a wide margin 
allowed for errors in these estimates. Since our arguments point to potentially different values 
of interest, in what follows we will allow for signihcant ranges of these Lovelock parameters. 

In order to simplify our analysis of the inhuence of the various Lovelock invariants 
it is perhaps best to examine them one at a time; we hrst set /? = 7 = 0 and consider the 
case of a 7 ^ 0. This is the most familiar possibility corresponding to Gauss-Bonnet gravity 
discussed in the literature [23 121 ] • In order to elucidate the effects of a 7 ^ 0 for any hxed 
value of n and Mbh it is useful to scale out the results by the predictions of the usual EH 
action as can be done using the expressions above. As we will see the various quantities of 
interest will have a reasonably simple dependence on a for hxed values of m = Mbh / M^. For 
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Figure 2: Same as in the previous figure but now for the BH temperature (top) and en- 
tropy(bottom). 
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Figure 3: Same as in the previous figure but now for the BH heat capacity(top) and free 
energy (bottom). 
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example, in the specific case of the Schwarzschild radius we will examine the ratio R/Ro to 
which we now turn and which is displayed in Fig.^ Here we see that the effect of a non-zero 
a on i? is rather mild for all values of n unless a is larger than ~ 10“^. Interestingly, on 
the right hand side of the plot we see that the cases with larger values of n appear to be 
significantly less influenced by a non-zero a than do those with smaller n. These variations 
in R can lead to significant modifications to the anticipated BH cross section. For example, 
if a ~ 1 and n = 2(6) the BH production cross section, ~ would be reduced by more 
than a factor of 100(5) due to the presence of this new Lovelock in the action in comparison 
to EH expectations. 



Figure 4: Same as Fig.l but now as a function of (3 with a = 7 = 0 and n = 2 ^ n = 3 
since £3 vanishes when n = 2. 


Fig. El shows the corresponding results for the temperature and entropy ratios, T/Tq 
and S'/S'o, as functions of a. The results shown in these two plots appear to be reasonably 
(anti-)correlated as a result of of the entropy definition, Eq.(15). In the top panel we see that 
T/Tq initially decreases as a increases, then goes through a minimum near a ~ 0.03 — 0.08, 
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then increases with a rapidity which is highly n dependent. We see that as n increases 
the steepness of this temperature rise rapidly decreases. Again, larger n appear to be less 
sensitive to appreciable values of a. The reverse behavior is seen in the case of the entropy 
variation with both a and n. Note that in both cases the ratios never differ from unity by 
more than factors of a few unless a becomes quite large. Thus, qualitatively, these BH are 
not too dissimilar from those of EH except for extreme choices of the parameter a. 



Figure 5: BH mass dependence of the Schwarzschild radius for the case a = 
7 = 0. Solid(dashed,dotted) curves correspond to n = 3(4,6), respectively. The 
red(green,blue,magenta,cyan) curves correspond to j3 = 0.00003(0.0001, 0.0003, 0.001, 0.003), 
respectively. Note here that R is in units of M~^. 


In Fig. 01 we see the behavior of both the heat capacity and free energy ratios with 
variations in a. The values of these two quantities are reasonably simple and qualitatively 
similar to both T and S. F/Fq has an easily understandable behavior with increasing a since 
F = m — TS and given the behavior above of the ratios T/Tq and S/Sq. It is interesting 
to note that at large a the ratio F/Fq asymptotes to a set of values which depend on n 
and not m. As in the EH case, dF/dm can be shown to be always positive and the heat 
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capacity for BH is seen to be well behaved and is always negative as in the EH case. Thus BH 
with only non-zero a will Hawking radiate, evaporating completely or until a Planck-scale 
remnant stage is reached. For very large values of a ~ 1, the mass loss through evaporation 
may be somewhat slowed in comparison to the usual EH case since the BH heat capacity, 
dm/dT, is so small. From Figs. 1-3, however, it is generally apparent that when n > 2 
and only a is non-zero the typical BH will behave in a manner which is qualitatively similar 
to (though quantitatively different from) those that arise solely from the EH action. They 

are all unstable and will have relatively short lifetimes typical of more familiar hadronic 
processes. 

The behavior of BH properties with non-zero (3 can be signihcantly different than 
what happens in the case of non-zero a as we will now see. (For simplicity we will assume 
a = 7 = 0 while studying non-zero (3; we note, however, that the behavior of non-zero 7 is 
qualitatively similar to non-zero /3.) We can see this as soon as we examine the (3 dependence 
of the ratio R/Rq as is shown in Ref. |31 For small values of (3 we see that the ratio R/Rq 
decreases away from unity with smaller values of n most seriously affected. This is quite 
similar to what we saw happen in the the case of non-zero a. Clearly either large a or large 
(3 can lead to a signihcant suppression of the BH production cross section; the same will 
be true of non-zero 7 . Here, however, we see something new for the case n = 3: the ratio 
R/Rq —>■ 0 for hnite (3. Turning this around, this means that for n = 3 a non-zero (3 implies 
that an horizon will not form unless m exceeds a certain critical value, he., a BH will not 
form unless a certain mass threshold is reached. (Such a possible behavior has already been 
observed within the general context of Lovelock extended actions^!?! I24ji. This threshold 
effect can be seen analytically from Eq.(lO) by setting n = 3 with a = 7 = 0 which yields 
the simple relation 

m = c{x‘^ + 24/?), (20) 
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Figure 6 : Same as Fig.2 but now as a function of (3 with a = 7 = 0 and with n = 2 


n = 3. 
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where here c = 57r^/2 and which is trivial to invert analytically to 


X = 



(21) 


Here we see that nnless m > there is no horizon. Note that for 13 = 5 ■ 10“^, a typical 

valne, this means that nnless m > 0.93 no BH are formed. It is important to observe that 
a similar resnlt is not obtainable when n > 3; the case n = 3 is therefore special when 
[3^0. Note that for hxed (3 simnltaneonsly tnrning on a non-zero a will lead to this same 
qnalitative behavior with the same minimnm valne of m. Perhaps even more interesting is 
the fact that something almost identical happens in the case of non-zero 7 when n = 5. 

This horizon threshold effect for non-zero valnes of (3 can perhaps be best seen in 
Fig. El where we show i? as a fnnction of m = Mbh/M^. Here it is qnite evident that 
there is a hnite horizon threshold for the case n = 3 while none exists for other n valnes, 
he., i? —>■ 0 when m —0. For n > 3 this is easily nnderstood from Eq.lO as there we 
see that m = 0 when x = 0. It is interesting to note that this threshold effect is not too 
dissimilar from what may happen for a 4-d BH when a renormalization gronp rnnning of 
Newton’s constant is employedin order to approximate leading qnantnm corrections; 
snch a threshold scenario can also be seen to occnr in theories with a minimnm length |2Z] 
and in loop qnantnm gravityj^H], which may also be thonght of as a qnantnm effect. It is 
interesting to note that this threshold phenomena occnrs in all these models where one tries 
to incorporate qnantnm corrections in some way; it would be interesting to learn whether or 
not this is a general feature of all such models. 

We might expect that this unusual threshold-type behavior will have a signihcant 
influence on thermodynamic quantities. Fig. El shows the dependence of T/Tq and S/Sq as 
functions of {3 for different values of n and m. For n > 3 this behavior is not so qualitatively 
different than that seen in the case of a non-zero a which is consistent with what was observed 
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Figure 7: Same as Fig. 5 but now for the BH temperature and entropy with u = 2 —> n = 3 
as above. Note that here the temperature is in units of M*. 
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Figure 8 : Same as Fig.3 but now as a function of (3 with a = 7 = 0. 
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in Figs. 0 and El However, in Fig. El we see that for n = 3 both the temperature and entropy 
can vanish. These zeros occur precisely at the values at which R ^ 0 for the n = 3 case that 
we observed above: T, S and R all vanish at the same parameter space points. These zeros 
in both the temperature and entropy for hnite m are also observed for the case of n = 3 
in Fig. El but we also see something new for other values of n. While dS/dm was always 
positive, it is possible for dT/dm to change sign in the mass range of interest. In the EH 
case, as well as for the examples of non-zero a examined above, dT/dm was always seen to 
be negative. In many cases for /? 7 ^ 0, T{m) is hrst observed to rise, reach a maximum, then 
decrease as m grows larger. 

To understand this more fully we turn to Figs. |H1 and E which show the BH heat 
capacity and free energy both of which can now have rather complex behavior. The least 
unusual of the two is the ratio F/Fq, which apart from the case of n = 3, is qualitatively 
similar to the that seen above for the case of a non-zero a. However, we still see that F can 
become negative for some values of jS. On the otherhand, the specihc heat/heat capacity is 
observed to have as many as two sign reversals associated with singularities in the parameter 
range of interest. In most cases C/Cq is positive for both large and small (3 with the reversals 
taking place at intermediate values. After some algebraic manipulations, it can be shown 
that for hxed values of m, C{(3)/Cq generally has two simple poles for real positive (3 values. 
Such behavior is not unheard of for BH arising from Lovelock extended actions piU and 
may even occur in 4-d|2n] for Kerr-Neumann BH. For the n = 3 case, however, since the 
range of R and m is restricted there can be only one singularity. 

If we instead examine C{m), as shown in Fig. El we see that there is at most one 
singularity corresponding to a simple pole in the physical region. For small m we see that 
C > 0 and then flips sign as m increases above the singularity. This signals the existence of 
a kind of phase transition occurring at this singularity. F{m) is also unusually behaved: in 
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^bh/^* 


Figure 9: Same as the Fig. 7 but now for the BH heat capacity and free energy. Recall that in 
our notation F is units of M* while the heat capacity is dimensionless. Curves terminating 
(and not changing sign) in the middle of the plot correspond to locations where i? ^ 0. 
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Figure 10: Threshold behavior of the BH cross section at the LHC(top) and ILC(bottom) for 
n = 3. The red curve assumes the absence of Lovelock terms, while the green(blue,magenta) 
one corresponds to = 0.0005 with a = 0(0.01,0.02). In the LHC case a luminosity (bin 
width) of 100 fb~^ (100 GeV) has been assumed. In the ILC case a scaled cross section is 
presented as discussed in the text; the usual step-function turn on in the absence of Lovelock 
terms has been assumed. For the ILC -y/i = Mbh- 
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general it is positive for both small and large m but goes through a negative minimum in 
between. Not unsurprisingly the location of the minimum in the free energy is the same as 
that of the singularity in the specihc heat/heat capacity. 

A parameter regime with positive heat capacity leads to interesting properties for BH. 
With only the EH action and thus a negative C, a BH loses mass via Hawking radiation; as 
it does so it get hotter thus radiating more quickly. This run-away process stops only when 
evaporation is complete or a non-classical Planck scale remnant is reached. Here, with /3 7 ^ 0, 
it is possible that the BH will end up in a mass regime where C* > 0 in which case it will get 
cooler as it gets lighter. Furthermore, there is a maximum temperature, given by the location 
of the pole in the heat capacity, for the BH. One might expect that such BH will live much 
longer than naively expected based on the EH action or even the combined EH+/I 2 action. 
This is something completely new that happens when higher order Lovelock invariants are 
present. In some cases this may imply that the BH will become classically stable which is 
what happens in the case of n = 3. To see this analytically, it is interesting to calculate the 
BH lifetime (in a vacuum!) directly assuming the semi-classical approach is viable. Recall 
that the amount of energy emitted by a black body per unit time is proportional to both 
its area and the fourth power of its temperature; using the results presented in Ref. j25j for 
n = 3 with a = 7 = 0 as an example we obtain 


dm 

dt 




(m — 

(m -|- 1207r3/3)^ 


( 22 ) 


Integrating from some arbitrary initial mass, mo, down to the threshold value, bOvr^/? (where 
T = 0), we obtain an inhnite lifetime for all non-zero values of /?. By ‘inhnite lifetime’ we 
mean that it will take an inhnitely long time for the mass of the BH to reach the threshold 
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value. These mass thresholds may be identihable with the anticipated Planck-scale remnant. 
We note that if /? = 0 a hnite lifetime would be obtained. The addition of a non-zero value 
of a will not change this inhnite lifetime result in a qualitative way. Thus any BH with a 
mass above the threshold value will decay by Hawking radiation quite rapidly losing most 
of its ‘excess’ mass until its mass approaches the threshold. Interestingly, we see from the 
hgures that as the BH loses mass and approaches the threshold mass value from above the 
free energy approaches a minimum which is another indicator of classical stability. Of course, 
taking (3^0 one recovers the hnite lifetime that is usually expected on the basis of the 
EH action. This apparent BH stability is of course not unique to these particular parameter 
values; other BH do not lead to hnite mass hnal states but may take an inhnite time to 
evaporate completely. Furthermore, for n = 5, we hnd that a similar situation is seen to 
occur when 7 is non-zero for arbitrary positive values of a, (3. Thus the so-called BH remnant 
in these cases can be uniquely identihed with a BH with a mass very close to the threshold 
value. 

At this point it is important that the reader be reminded that even though we are 
trying to model quantum/stringy corrections via the Lovelock invariants, these are still semi- 
classical predictions which are subject to true quantum corrections that may modify them 
signihcantly. Though classically stable, these BH may have hnite lifetimes due to these 
quantum corrections. It is possible, however, that these lifetimes will remain quite long. 

4 A Toy Model of BH Collider Cross Section Thresh¬ 
olds 

As discussed above, the usual picture of collider BH production is that there is a 6 ^ function- 
type threshold at M* which is certainly unphysical. Clearly quantum and/or stringy correc¬ 
tions will turn this step function behavior into something smooth which might be modeled 
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Figure 11: Close up of the BH production threshold at the ILC for n = 3(top) and n = 
5(bottom); in these plots we should remember that ^/s = Mbh- The blue line in the naive 
6 '-function that is usually assumed in the absence of Lovelock terms. For the case n = 3 with 
/3 = 0.0005 the curves correspond to a = 0, 0.002, 0.005, 0.01 and 0.02 from top to bottom, 
respectively. For the case n = 5 all curves have 7 ~ 1.14 ■ 10“^ and, from top to bottom, 
have a, /3 = 0, 0(0.003, 0; 0, 0.00003; 0.01, 0; 0, 0.0001), respectively. 
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via additional Lovelock terms in the action. If TeV scale extra dimensions exist and BH are 
produced at future colliders, the threshold region, where these effects are most relevant, will 
be important to understand in detail since it will provide a handle on the true underlying 
theory of quantum gravity. In the scenario we propose below, this corresponds to extracting 
the values of the Lovelock parameters by examining the nature of the cross section threshold. 
This is akin to using the nature of the threshold for new particle production to determine 
the particle’s spin and other quantum numbers. 

The analysis in the previous section suggests a toy model that yields a more properly 
behaved threshold for BH production, at least in the cases n = 3(5), when /^(y) is non-zero. 
There we saw that, e.g., a non-zero (3 for n = 3 implies the existence of a minimum mass 
below which a BH simply will not form. We also saw that this mechanism was absent for even 
values of n. For reasonable values of the Lovelock parameters what do these BH thresholds 
look like? Sticking with the case n = 3 for now, let us take (3 = 0.0005, in the center of 
the range discussed above with a set to zero for the moment for simplicity. Note that since 
n = 3, 7 is automatically zero. The top panel in Fig. shows what this would imply 
for the LHC BH (very large) production rate assuming that M* = 1.5 TeV for purposes 
of demonstration. First we notice the red histogram which is the usual ^-function turn-on 
cross section prediction. The green histogram shows the same cross section but now with 
(3 = 0.0005 and no 6*-function. There is clearly a threshold near 1.5 TeV but the turn-on 
is no longer inhnitely sharp; for large BH masses we see that this prediction asymptotes to 
that of the naive ^-function histogram. This result we may have expected since for large 
values of Mbh we should recover the predictions associated with the EH action. Turning on 
non-zero but typical values of a ~ 0.01 does not alter the threshold location but does modify 
the cross section shape near the threshold. Apparently the larger the value of a the softer 
the slope of the cross section near threshold. In fact a short analysis shows that the slope of 
the cross section goes as ~ 1/a precisely at the threshold. Of course, radiative corrections 
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will smear out the fine details of these rather sharp threshold shapes. For non-zero a the 
cross section is seen to still asymptote to the naive value in the large Mbh limit as it should. 

Although a detailed simulation is far beyond the scope of the present work it is clear 
from the top panel in Fig.Elthat, while the LHC will observe the BH production threshold, 
an extraction of the Lovelock parameter values with any precision will be difficult. The 
bottom panel of this figure shows what may be seen at the ILC; here a set of scaled cross 
section results are presented so that the reader can move the threshold as they prefer. For 
example, if M^, = 1 TeV, crM^ = 1 corresponds to a cross section of ~ 389.4 pb; there 
will be no shortage of statistics. As in the case of the LHC the red curve is the naive step 
function turn-on while the other curves show the effect of the non-zero Lovelock parameters. 
To get a better idea of the ILC threshold region we zoom in using the top panel of Fig. 

In this idealized picture we see that the step function curve is very far away from any of 
those produced by the Lovelock terms. The Lovelock cross section predictions are themselves 
rather widely separated from each other and have quite distinct shapes as they approach 
the threshold. Again, we recall that if M* = 1 TeV, a single vertical tickmark on this 
plot corresponds to a cross section change of ~ 8 pb, he., ~ 10® events at the ILC in one 
year. Of course in the real world this threshold region will be made more complex by both 
electroweak and gravitational radiative corrections that may smear out the threshold to some 
extent. Just how significant such effects will be is not yet known even within the very limited 
scope of the present toy model. 

For the case of n = 5 a non-zero 7 also induces a BH mass threshold and the bottom 
panel shows the corresponding cross sections assuming 7 ~ 1.1 ■ 10 “®, a value in the antici¬ 
pated region discussed above. Note that the slope of the cross section in the threshold region 
is generally increased as n increases. Turning on different values of a and (3 again leads to 
detailed modifications of the threshold cross section shape. When a = /? = 0, the slope of 
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the cross section at threshold is inhnite; a short analysis shows that in this case the slope 
behaves as ~ 1//3 near threshold and is independent of a. The various curves are seen to 
be sufficiently different as to be likely distinguishable at the ILC even after radiative effects 
are included. Thus we see that the cases P ^ 0 and 7 7^ 0 are qualitatively similar but are 
clearly separable once measurements with any precision become possible. 

At least for the case of n odd we see that the Lovelock invariants can lead to a 
reasonably smooth turn on of BH production in the mass range near M*. The precise nature 
of this turn in is highly correlated with the various parameter values and the number of extra 
dimensions. 


5 Conclusions 

In this paper we have examined how the production of TeV-scale black holes at colliders is 
influenced by the existence of higher curvature Lovelock invariants in the action of models 
with large extra dimensions. This possibility is anticipated on rather general grounds when 
we try to embed such models in a larger theoretical structure, e.g., strings, if we also wish to 
maintain unitarity and freedom from ghosts. In the case of the ADD model, we argued that 
the addition of the Lovelock invariants to the Einstein-Hilbert action would leave the more 
conventional missing-energy and spin-2 contact interaction signatures unaltered but would 
modify the anticipated properties of TeV scale BH since the ADD geometry is essentially 
flat. These modifications in BH properties provide the only signature for the existence of 
higher curvature terms in this context. Since the total number of dimensions D = A + n 
in string-motivated ADD is < 10 and the properties of the Lovelock invariants are highly 
restrictive, it was easily shown that at most three such terms, ^^ 2 , 3 , 4 , can be added to the 
usual Einstein-Hilbert action. The strengths of these new interactions are described by 
a set of dimensionless parameters and, given the expectation that such terms result from 
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some form of perturbative expansion, we derived some rough estimates for their anticipated 
magnitudes. 

The next step in our analysis was to obtain the general expressions for the Schwarzschild 
radius, Hawking temperature, entropy, heat capacity/specihc heat and free energy relevant 
for the D-dimensional TeV-scale BH in Lovelock gravity models with our specihc generalized 
action. Using these results we then performed a detailed numerical investigation to probe 
the sensitivity of the BH production cross section as well as these thermodynamic quanti¬ 
ties to variations in the values of the Lovelock parameters. Since the number of potentially 
contributing invariants is n-dependent, we examined the possibility of non-zero values for 
coefficients (a and {3) of both £ 2,3 one at a time; while the impact of £2 and £3 were found to 
be quite distinct, £4 leads to results which are quite similar to £3 and so was not examined 
in detail. 

For non-zero values of a, up to factors or order unity, the BH thermodynamical 
properties are overall qualitatively similar to those found in the pure EH case except when 
a get very large. The BH Schwarzschild radius is, however, signihcantly smaller which 
suppresses the collider BH production cross section. When (3 is non-zero, both quantitative 
as well as qualitative variations in the BH properties can be anticipated especially for the case 
of n = 3; these variations from EH expectations can be somewhat exotic. We showed that 
they include sign changes in the heat capacity/specihc heat, indicative of phase transitions, 
negative free energies and BH mass thresholds, he., mass values below which BH cannot 
form. In particular we observed that BH can now become thermodynamically stable, he., 
it would take an inhnite time for them to Hawking radiate down to their threshold masses. 
This is correlated with BH now possibly having positive heat capacities so that they cool as 
they lose mass thus slowing the mass loss process itself. 

We then demonstrated that for n = 3, 5 the minimum BH masses induced by the 
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Lovelock invariants could be used as a toy model for the BH production cross section thresh¬ 
old. Instead of a step function turn-on, the cross section near threshold was shown to be 
highly sensitive to the values of a, and 7 yet would still asymptote to that expected in the 
pure EH case. These differences would surely show up at the LHC and can be probed in 
detail at the ILC due to the large cross section and the enormous statistics available. It may 
be possible to use these BH mass thresholds to extract the values of the Lovelock parameters 
themselves. 

TeV scale BH can provide a unique tool to probe the underlying theory of gravity. 
Hopefully they will be seen at future colliders. 
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